home *** CD-ROM | disk | FTP | other *** search
/ Java Programmer's Toolkit / Java Programmer's Toolkit.iso / gs3.53 / docie.ps < prev    next >
Text File  |  1996-01-10  |  7KB  |  220 lines
  1. %    Copyright (C) 1995 Aladdin Enterprises.  All rights reserved.
  2. % 
  3. % This file is part of Aladdin Ghostscript.
  4. % 
  5. % Aladdin Ghostscript is distributed with NO WARRANTY OF ANY KIND.  No author
  6. % or distributor accepts any responsibility for the consequences of using it,
  7. % or for whether it serves any particular purpose or works at all, unless he
  8. % or she says so in writing.  Refer to the Aladdin Ghostscript Free Public
  9. % License (the "License") for full details.
  10. % 
  11. % Every copy of Aladdin Ghostscript must include a copy of the License,
  12. % normally in a plain ASCII text file named PUBLIC.  The License grants you
  13. % the right to copy, modify and redistribute Aladdin Ghostscript, but only
  14. % under certain conditions described in the License.  Among other things, the
  15. % License requires that the copyright notice and this notice be preserved on
  16. % all copies.
  17.  
  18. % docie.ps
  19. % Emulate CIE algorithms in PostScript.
  20.  
  21. % ---------------- Auxiliary procedures ---------------- %
  22.  
  23. /r1default [0 1] def
  24. /r3default [0 1 0 1 0 1] def
  25.  
  26. /apply3        % <u> <v> <w> [<pu> <pv> <pw>] apply3 <u'> <v'> <w'>
  27.  {  { 4 -1 roll exch exec } forall
  28.  } bind def
  29.  
  30. /restrict    % <u> <min> <max> restrict <u'>
  31.  { 3 1 roll .max .min
  32.  } bind def
  33.  
  34. /restrict3    % <u> <v> <w> [<minu> ... <maxw>] restrict3 <u'> <v'> <w'>
  35.  { aload pop
  36.    7 -1 roll 3 1 roll restrict 7 1 roll
  37.    5 -1 roll 3 1 roll restrict 5 1 roll
  38.    restrict 3 1 roll
  39.  } bind def
  40.  
  41. /rescale    % <u> <min> <max> rescale <u'>
  42.  { 1 index sub 3 1 roll sub exch div 0 .max 1 .min
  43.  } bind def
  44.  
  45. /rescale3    % <u> <v> <w> [<minu> ... <maxw>] rescale3 <u'> <v'> <w'>
  46.  { aload pop
  47.    7 -1 roll 3 1 roll rescale 7 1 roll
  48.    5 -1 roll 3 1 roll rescale 5 1 roll
  49.    rescale 3 1 roll
  50.  } bind def
  51.  
  52. /mmult3        % <u> <v> <w> [<uu> <uv> ... <wv> <ww>] mmult3
  53.         %   <u'> <v'> <w'>
  54.  { 4 -1 roll dup dup 6 -1 roll dup dup 8 -1 roll dup dup
  55.    10 -1 roll { 10 -1 roll mul } forall
  56.         % Stack: u1 v1 w1 u2 v2 w2 u3 v3 w3
  57.    4 -1 roll add 6 -1 roll add
  58.         % Stack: u1 v1 u2 v2 u3 v3 w'
  59.    7 1 roll 3 -1 roll add 4 -1 roll add
  60.         % Stack: w' u1 u2 u3 v'
  61.    5 1 roll add add 3 1 roll
  62.  } bind def
  63.  
  64. /minvert3    % [<uu> <uv> ... <wv> <ww>] minvert3
  65.         %   [<uu'> <uv'> ... <wv'> <ww'>]
  66.  { 16 dict begin
  67.    aload pop { I H G F E D C B A } { exch def } forall
  68.    /coa E I mul F H mul sub def
  69.    /cob F G mul D I mul sub def
  70.    /coc D H mul E G mul sub def
  71.    /det A coa mul B cob mul add C coc mul add def
  72.    [ coa det div
  73.      C H mul B I mul sub det div
  74.      B F mul C E mul sub det div
  75.      cob det div
  76.      A I mul C G mul sub det div
  77.      C D mul A F mul sub det div
  78.      coc det div
  79.      B G mul A H mul sub det div
  80.      A E mul B D mul sub det div
  81.    ]
  82.    end
  83.  } bind def
  84.  
  85. /print1
  86.  { print dup ==
  87.  } bind def
  88.  
  89. /print3
  90.  { print 3 array astore dup == aload pop
  91.  } bind def
  92.  
  93. % ---------------- Mapping to XYZ ---------------- %
  94.  
  95. /csmap        % <csdict> <l> <m> <n> csmap <csdict> <x> <y> <z>
  96.  { 3 index /RangeLMN .knownget not { r3default } if restrict3
  97. DEBUG { (After RangeLMN Decode: ) print3 } if
  98.    3 index /DecodeLMN .knownget { apply3 } if
  99. DEBUG { (After DecodeLMN Decode: ) print3 } if
  100.    3 index /MatrixLMN .knownget { mmult3 } if
  101. DEBUG { (After MatrixLMN Decode: ) print3 } if
  102.  } bind def
  103.  
  104. /csciea        % <csdict> <a> csciea <csdict> <x> <y> <z>
  105.  { 1 index /RangeA .knownget not { r1default } if restrict
  106. DEBUG { (After RangeA Decode: ) print1 } if
  107.    1 index /DecodeA .knownget { exec } if
  108. DEBUG { (After DecodeA Decode: ) print1 } if
  109.    1 index /MatrixA .knownget
  110.     { { 1 index mul exch } forall pop }
  111.     { dup dup }
  112.    ifelse
  113. DEBUG { (After MatrixA Decode: ) print3 } if
  114.    csmap
  115.  } bind def
  116.  
  117. /cscieabc    % <csdict> <a> <b> <c> cscieabc <csdict> <x> <y> <z>
  118.  { 3 index /RangeABC .knownget not { r3default } if restrict3
  119. DEBUG { (After RangeABC Decode: ) print3 } if
  120.    3 index /DecodeABC .knownget { apply3 } if
  121. DEBUG { (After DecodeABC Decode: ) print3 } if
  122.    3 index /MatrixABC .knownget { mmult3 } if
  123. DEBUG { (After MatrixABC Decode: ) print3 } if
  124.    csmap
  125.  } bind def
  126.  
  127. % ---------------- Rendering from XYZ ---------------- %
  128.  
  129. /lookup3    % <rtable> <a[0..1]> <b[0..1]> <c[0..1]> lookup3
  130.         %   <rtable> <bytes>
  131.  { 3 -1 roll 3 index 0 get 1 sub mul
  132.    3 -1 roll 3 index 1 get 1 sub mul
  133.    3 -1 roll 3 index 2 get 1 sub mul
  134.         % Stack: rtable ia ib ic
  135. DEBUG { (RenderTable indices: ) print3 mark 5 1 roll } if
  136.    3 -1 roll round cvi 3 index 3 get exch get
  137.         % Stack: rtable ib ic string
  138.    3 -1 roll round cvi 3 index 2 get mul
  139.         % Stack: rtable ic string ib*nc
  140.    3 -1 roll round cvi add 2 index 4 get mul
  141.         % Stack: rtable string index
  142.    2 index 4 get getinterval
  143.         % Stack: rtable bytes
  144. DEBUG { (RenderTable values: ) print (<) print (%stdout) (w) file 1 index writehexstring (>) = } if
  145.  } bind def
  146.  
  147. /bpdefault [0 0 0] def
  148.  
  149. /crmap        % <csdict> <crdict> <x> <y> <z> crmap <v1> ...
  150.  {
  151. DEBUG { (CIE XYZ = ) print3 } if
  152.    3 index /MatrixPQR .knownget { mmult3 } if
  153. DEBUG { (After MatrixPQR: ) print3 } if
  154.    4 index /WhitePoint get
  155.    5 index /BlackPoint .knownget not { bpdefault } if
  156.    5 index /WhitePoint get
  157.    6 index /BlackPoint .knownget not { bpdefault } if
  158.    4
  159.      { 4 -1 roll aload pop
  160.             % Stack: csdict crdict x y z pt pt pt px py pz
  161.        3 copy 12 index /MatrixPQR .knownget { mmult3 } if 6 array astore
  162.      }
  163.    repeat
  164.             % Stack: csdict crdict x y z wps+ bps+ wpd+ bpd+
  165.    9 -1 roll pop        % get rid of csdict
  166.    7 4 roll
  167.    7 index /TransformPQR get
  168.     {            % Stack: crdict wps+ bps+ wpd+ bpd+ u v w proc
  169.       8 copy exch pop exch pop
  170.       exec exch pop 4 -1 roll pop
  171.     }
  172.    forall
  173.    7 3 roll pop pop pop pop    % get rid of White/BlackPoints
  174. DEBUG { (After TransformPQR: ) print3 } if
  175.    3 index /MatrixPQR .knownget { minvert3 mmult3 } if
  176. DEBUG { (After MatrixPQR': ) print3 } if
  177.    3 index /MatrixLMN .knownget { mmult3 } if
  178. DEBUG { (After MatrixLMN Encode: ) print3 } if
  179.    3 index /EncodeLMN .knownget { apply3 } if
  180. DEBUG { (After EncodeLMN Encode: ) print3 } if
  181.    3 index /RangeLMN .knownget not { r3default } if restrict3
  182. DEBUG { (After RangeLMN Encode: ) print3 } if
  183.    3 index /MatrixABC .knownget { mmult3 } if
  184. DEBUG { (After MatrixABC Encode: ) print3 } if
  185.    3 index /EncodeABC .knownget { apply3 } if
  186. DEBUG { (After EncodeABC Encode: ) print3 } if
  187.    3 index /RangeABC .knownget not { r3default } if
  188.    5 -1 roll /RenderTable .knownget
  189.     {        % Stack: u v w ranges rtable
  190.       5 1 roll rescale3
  191. DEBUG { (Rescaled ABC: ) print3 } if
  192.         % Stack: rtable a b c
  193.       lookup3
  194.         % Stack: rtable bytes
  195.       0 1 3 index 4 get 1 sub
  196.        {    % Stack: values rtable bytes c
  197.      2 copy get 255 div
  198.         % Stack: values rtable bytes c v
  199.      3 index 3 -1 roll 5 add get exec 3 1 roll
  200.        }
  201.       for pop pop
  202. DEBUG { (After RenderTableT: ) print ] dup == aload pop } if
  203.     }
  204.     { restrict3
  205. DEBUG { (After RangeABC Encode: ) print3 } if
  206.     }
  207.    ifelse
  208.  } bind def
  209.  
  210. % ---------------- Top level control ---------------- %
  211.  
  212. /mapdict mark
  213.   /CIEBasedA { 1 get exch csciea currentcolorrendering 4 1 roll crmap } bind
  214.   /CIEBasedABC { 1 get 4 1 roll cscieabc currentcolorrendering 4 1 roll crmap } bind
  215. .dicttomark def
  216. /mapcie        % <a> mapcie <v1> ...
  217.         % <a> <b> <c> mapcie <v1> ...
  218.  { currentcolorspace dup 0 get //mapdict exch get exec
  219.  } bind def
  220.